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Abstract. Various satellite-borne missions are being planned to measure the 
polarization of a large number of gamma-ray bursts (GRBs). We show that the 
polarization pattern resulting from the current models of GRB emission can be 
drastically modified by the existence of very light axion-like particles (ALPs), which 
are predicted by many extensions of the Standard Model of particle physics. Basically, 
the propagation of photons emitted by a GRB through cosmic magnetic fields with 
a domain-like structure induces photon-ALP mixing, which is expected to produce 
a strong modification of the initial photon polarization. Because of the random 
orientation of the magnetic field in each domain, this effect strongly depends on 
the orientation of the line of sight. As a consequence, photon-ALP conversion 
considerably broadens the initial polarization distribution. Searching for such a 
peculiar feature through future high-statistics polarimetric measurements therefore 
offers a new opportunity to discover very light ALPs. 
Keywords: axions, gamma-ray burst polarization. 
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1. Introduction 

A generic feature of many extensions of the Standard Model is the prediction of 
axion-like particles (ALPs), namely very light spin-zero bosons characterized by a two- 
photon coupling 077. Besides occurring in four-dimensional models [1, 2, 3, 4], ALPs 
naturally arise within compactified Kaluza-Klein theories [5] as well as in superstring 
theories [6, 7]. As the name itself suggests, ALPs are a generalization of the axion, 
the pseudo-Goldstone boson associated with the Peccei-Quinn symmetry proposed as a 
natural solution to the strong CP problem [8, 9, 10, 11, 12]. Specifically, while the axion 
is characterized by a strict relationship between its mass and 077 coupling constant, 
these two parameters are to be regarded as independent for ALPs [13]. Depending on 
the actual values of their mass and 077 coupling constant, ALPs can play an important 
role in cosmology, either as cold dark matter particles responsible for the structure 
formation in the Universe [14] or as quintessential dark energy [15] which presumably 
triggers the present accelerated cosmic expansion (see [16] for a recent review on ALPs). 

Owing to the 077 coupling, in an external electromagnetic field the phenomenon 
of photon- ALP mixing takes placef. More specifically, two very interesting effects 
arise in such a situation. One is photon-ALP conversion (oscillation) [17, 18, 19], 
which is exploited by the AD MX experiment to search for ALP dark matter [20], 
by CAST to search for solar axions [21, 22, 23] and by the regeneration laser 
experiments [24, 25, 26, 27, 28, 29, 30]. The other effect consists in the change of 
the polarization state of photons traveling in a magnetic field. In particular, an initially 
linearly polarized photon beam propagating in a transverse magnetic field acquires an 
elliptical polarization with the major axis rotated with respect to the direction of the 
initial polarization [17, 31, 32]. A claim for a positive observational evidence of such an 
effect by the PVLAS collaboration [33] employing a laser beam has subsequently been 
withdrawn [34]. 

Laboratory experiments devised to search for an ALP-induced photon beam 
polarization suffer from the intrinsic limitation that a very short baseline is available. 
Remarkably enough, astrophysical observations of distant X-ray and 7-ray sources offer 
new opportunities to look for ALP effects in polarization measurements (see [35] for 
a recent comprehensive study). Actually, several satellite-borne missions are currently 
under consideration to perform the challenging measurement of the polarization state 
of photons emitted by distant astrophysical X-ray and 7-ray sources like the gamma-ray 
bursts (GRBs) in the keV-MeV energy range. These new missions have a wide field 
of view and a broad energy band, thereby guaranteeing polarimetric measurements of 
GRBs to be performed with high statistics. Such a circumstance evidently allows for a 
reliable determination of the statistical properties of the GRB polarization, which has 
recently been proposed as a crucial tool to discriminate among different models for the 
GRB emission [36]. 

I We stress that the external field is necessary in order to compensate for the photon-ALP spin 
mismatch. 
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We show that strong modifications in the polarization pattern of distant X-ray and 
7-ray sources can be produced by ALPs with parameters lying in experimentally allowed 
ranges. Basically, the following situation is envisaged. Photons are emitted by the 
considered sources at cosmological distances and on their way to us they cross different 
magnetic field configurations, which ought to induce photon- ALP mixing. Consequently, 
the initial photon polarization gets changed and this change can eventually be detected. 

Manifestly, the role of the magnetic field is crucial in this respect and we have 
to contemplate all kinds of magnetic field configurations that the photon beam can 
experience, namely the magnetic field inside the source, the extragalactic magnetic 
field, the Galactic magnetic field and possibly an intracluster magnetic field if the 
beam crosses a cluster of galaxies. All these magnetic field components have a quite 
complicated and poorly known morphology, so that it has become customary to suppose 
that they possess a domain-like structure with varying coherence lengths. Therefore, 
we are actually dealing with a long baseline astrophysical setup to study ALP effects 
on the polarization of X-ray and 7-ray photons from distant sources. Our aim is indeed 
to investigate under which circumstances these effects are sizeable and detectable. For 
definiteness, our attention will be focussed on GRBs but our conclusions apply to any 
far-away astrophysical source of photons in the keV-MeV energy range, among which 
blazars play a very important role. We stress that at variance with a previous analysis 
of similar effects concerning radio sources [37] our picture does not require the actual 
presence of an ALP cosmological background in the Universe but merely demands that 
ALPs are produced by photon oscillations. 

The plan of the paper is as follows. In Section 2 we discuss the perspectives and 
the expectations for the measurement of GRB polarization. In Section 3 we review 
the mechanism of photon- ALP mixing in random magnetic fields which is relevant for 
our physical case and we identify the allowed ranges of the parameters that determine 
the effect in question. In Section 4 we discuss the impact of photon-ALP mixing on 
the polarization of GRBs for various intervening magnetic field configurations. We find 
that photon-ALP mixing in the magnetic field of GRBs is strongly suppressed. Still, 
during photon propagation from the source to us a sizeable mixing can occur. Due to 
the randomness of crossed magnetic fields, photon-ALP mixing depends strongly on the 
orientation of the line of sight. As a consequence, this effect is expected to produce 
a peculiar broad distribution in the polarization of GRBs with different projected 
position on the sky. Thus, the detection of such a feature through future high-statistics 
polarimetric measurements can become a new strategy to look for the existence of ALPs 
with very low mass. Finally, in Section 5 we offer a discussion of our results as well as 
of the possible further developments in this field. 

2. Gamma-ray burst polarization measurements 

GRBs are brief, intense flashes of 7-rays originating at cosmological distances and they 
are the most luminous objects in the Universe. Presently available lightcurves and 
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spectral observational information fail to provide a unique answer concerning their 
emission mechanism, and the polarization measurements of their X-ray and 7-ray 
emission have been recognized as a crucial mean to shed light on the inner structure of 
GRBs, including the geometry and the physical processes occurring close to the central 
engine. 

To date, GRB polarization measurements have been carried out systematically 
only in the optical band and in the afterglow phase. Sensitive observational 
techniques for X-ray and 7-ray polarimetry are being developed and in the next 
few years various polarimetric missions can finally enable to measure the X- 
ray and 7-ray polarization of the GRB emission. The missions presently under 
consideration include POET (Polarimeters for Energetic Transients) [38], PoGO 
(Polarimeter of Gamma-ray Observer) [39], POLAR [40], GEMS (Gravity and 
Extreme Magnetism) [41], XPOL [42], GRIPS (Gamma- Ray Burst Investigation via 
Polarimetry and Spectroscopy) [43] and NHXM (New Hard X-ray Mission) [44]. 

In particular, POET is the only one that incorporates a broadband capability 
for measuring the prompt GRB emission and employs two different polarimeters, both 
with a wide field of view: GRAPE (Gamma-Ray Polarimeter Experiment) operating 
between 60keV and 500 keV and LEP (Low Energy Polarimeter) operating between 
2 keV and 500 keV. Two smaller mission concepts based on the POET instruments are 
already scheduled for launch: GRAPE and GRBP (Gamma-Ray Burst Polarimeter, a 
smaller version of LEP). Both instruments offer the opportunity for a first glimpse at the 
polarization of prompt GRB emission and in fact POET provides sufficient sensitivity 
and sky coverage to detect up to 200 GRBs in a two-years mission. 

In general, given a random distribution of viewing angles, a statistical study of the 
polarization properties of a large sample of GRBs should allow to discriminate among 
different emission models and can provide a direct diagnostic tool for the magnetic 
field structure, the radiation mechanism and the geometric configuration of GRB jets. 
Toma et al. [36] display the predictions for the distribution of the amount of linear 
polarization in three different models: synchrotron emission with ordered magnetic 
fields (SO), synchrotron emission in random magnetic fields (SR) and Compton-drag 
model (CD). Specifically, the ratio N m /Nd of the number N m of GRBs for which the 
degree of polarization can be measured to the number N d of GRBs that are detected, 
and the distribution of the degree of linear polarization II^, can be used as criteria. It 
turns out that if N m /N d > 30% and Tl L clusters between 0.2 and 0.7, then the SO model 
will be favored. If instead N m /N d < 15%, then both the SR and the CD model will be 
preferred. Finally, if several events with U L > 0.8 are observed, then the CD model will 
instead be singled out. 

We will show that the presence of ALPs mixing with photons in astrophysical 
magnetic fields can drastically affect such expected statistical distributions for the linear 
polarization of GRBs. 
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3. Photon-ALP mixing in random magnetic fields 



3.1. Equations of motion 

The Lagrangian describing the photon-ALP system is 



a 7 



(1) 



The QED Lagrangian for photons is 



1 a 2 
4 + 90mf 



where F^ v = (E, B) is the electromagnetic field tensor, F t 



2 e Mi / po- 



(2) 

F pa is its dual, a 

is the fine-structure constant and m e is the electron mass. Natural Lorentz-Heaviside 
units with with h = c = k B = 1 are employed throughout. The second term on the r.h.s. 
of Eq. (2) is the Euler-Heisenberg-Weisskopf (HEW) effective Lagrangian [17], which 
accounts for the one-loop corrections to classical electrodynamics. The Lagrangian for 
the noninteracting ALP field a is 



C a = - d^a d u a --m 2 a 2 
2 ^2 



(3) 



where m is the ALP mass. A general feature of ALP models is the CP-conserving 
pseudo-scalar two-photon coupling, so that the photon-ALP interaction is represented 
by the following Lagrangian [17] 



'07 



~ gajF^F^a = g ai E • B a 



(4) 



where g ai is the photon-ALP coupling constant (which has the dimension of an inverse 

I 1/2 1/2 

energy). It is always assumed g ai Gp and m <^ G F , with G F ~ 250 GeV 
denoting the Fermi scale of weak interactions. 

We shall be concerned throughout with a monochromatic photon/ALP beam of 
energy E propagating along the z-direction in the presence of a magnetic field B. 
Clearly, the beam propagation is described by the second-order coupled generalized 
Klein-Gordon and Maxwell equations arising from the Lagrangian in Eq. (1). However, 
since we are interested in the regime E ^> m the short-wavelength approximation can 
be applied successfully and turns the beam propagation equation into the following 
Schrodinger-like one [17] 

A y (z) = , (5) 

' V a ( z ) J 

where A x (z) and A y (z) are the two photon linear polarization amplitudes along the x 
and y axis, respectively, a(z) denotes the ALP amplitude and Ai represents the photon- 
ALP mixing matrix. 



i^ + E + M 

, dz 
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For our futher purposes, it is more convenient to work with the polarization density 
matrix 

/ A x (z) \ 

( A x (z) A y (z) a(z) )* (6) 



(7) 



Ay(z) 

which obeys the Liouville-Von Neumann equation 

associated with Eq. (5). We denote by T(z, zq) the transfer function, namely the solution 
of Eq. (5) with initial condition T(z ,z ) = 1. Then any solution of Eq. (7) can be 
represented as 

p(z) = T(z, zq) p(z ) T\z, zq) . (8) 

The mixing matrix Ai takes a simpler form if we restrict our attention to the case 
of a photon beam propagating in a single magnetic domain, where the magnetic field 
B is supposed to be homogeneous. We denote by B T the transverse magnetic field, 
namely its component in the plane normal to the beam direction. We can choose the 
y-axis along B T so that B x vanishes. Under these simplifying assumptions, the mixing 
matrix can be written as [45] 



(9) 

(10) 
(11) 
(12) 





( A ± 








\ 







A,, 








V 


A a7 


A a 


/ 



whose elements are [17] 

A,, = A pl + 3.5A QED , 

A ± = A p i + 2 A Q ed , 
1 



4 a7 



A n = -— 



m 
2E 



1.52 x 10" 



-7.8 x 10" 



10- 11 GeV 

2 



Bn 



m 



10" 13 eV/ V10 2 keV 



E 



10~ 9 G 
i 



Mpc" 1 



with 



to 



A p i = ~ -1.1 x 10 



^QED 



2E 
aE 

457T 



E 



Bx 

B rr 



10 2 keV 
4.1 x 10~ 16 



n P 



10- 7 



cm 



Mpc" 1 , 



Mpc -1 



E 



10 2 keW V10" 9 G 



B 1 



Mpc" 1 



(13) 

(14) 
(15) 



where n e is the electron density in the medium, uj^ = 47ian e /m e is the associated plasma 
frequency, B CT = m 2 /e ~ 4.41 x 10 13 G is the critical magnetic field and e denotes the 
electron charge. 
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3.2. Input parameters 

Before proceeding further, we find it convenient to discuss the ranges of the parameters 
entering our scenario. The strength of the widespread, all-pervading magnetic field in 
the extragalactic medium has to meet the constraint B <2.8 x 1CT 7 (L/Mpc)~ 1 / 2 G - 
where L denotes its coherence length - which arises by scaling the original bound from 
the Faraday effect of distant radio sources [46, 47] to the now much better known baryon 
density measured by the Wilkinson Microwave Anisotropy Probe (WMAP) mission [48]. 
Its coherence length is expected to lie in the range 1 Mpc < L < lOMpc [49]. 
The mean diffuse intergalactic plasma density is bounded by n e <2.7 x 1CT 7 cm -3 , 
corresponding to the WMAP measurement of the baryon density [48]. Recent results 
from the CAST experiment yield a direct bound on the photon- ALP coupling constant 
g aj ^8.8 x KT 11 GeV" 1 for m^0.02 eV [23], slightly better than the long-standing 
globular-cluster limit [50]. In addition, for m<10~ 10 eV a more stringent limit arises 
from the absence of 7-rays from SN 1987A, giving g a ^^l x 1CT 11 GeV -1 [51, 52] even 
if with a large uncertainty§. 



3.3. Mixing in a single magnetic domain 



Inside a single magnetic domain B is homogeneous. Then it is straightforward to check 
that Ai^ in Eq. (9) can be brought into a diagonal form 

/ D 1 \ 



D = 



V 










D 2 






£3 / 

by the similarity transformation 
D = WM {0) W ] 



with 



W = 



( 1 \ 

cos 9 sin 9 
— sin 9 cos 9 



where the mixing angle 9 is given by 



9 = - arctan 
2 



2A 



£17 



All - A n 



and explicitly we get 

D 1 = A ± , 

(Aii+A a ) 



+ 



(A a -A||) 2 + 4A 



1/2 



(16) 



(17) 



(18) 



(19) 

(20) 
(21) 



§ These bounds on ALPs can be relaxed if they have a chameleontic nature [53] , in which case the most 
stringent constraint comes from the observed starlight polarization and reads g a7 < 10~ 9 GeV _1 [35]. 
However, we do not commit ourselves to chameleontic ALPs. 
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D* 



(A||+A ) 



1 



(A a -A||) 2 + 4A 



1/2 



(22) 
(23) 



2 2 

So, we find that the transfer matrix is presently 

T {0) (z,z ) = W^e w{z ~ Za) W . 

Finally, the density matrix can be computed by means of the general formula in Eq. (8) 
The probability that a photon initially polarized along the y axis (^22(0) = 1 
converts into an ALP after a distance d then reads [17] 



p(0) 

127 



p 33 (d) = sm 2 29 sin 2 ( 



where the oscillation wave number is 



(A a 



+ 4A 



0,7 



1/2 



(24) 



(25) 



Experience with problems similar to the one considered here shows that it proves useful 
to define a low critical energy [54] 



E T . = 



E\A n - A 



pil 



25 |m 2 



pi I 



10- y G\ /10- n GeV 



2 A, 



07 



(10" 13 eV) 2 V B T 



9a-y 



keV (26) 



along with a high critical energy [55] 
Eh 



~ 2 .1 x 10 18 



9a-y 



keV . 



(27) 



7 a B T 

It is easy to see that the oscillation wave number can be expressed in terms of El and 
Eh as 



2A 



0,7 



1 



sgn(m a - Wpij 



2^ 1/2 



(28) 



and from Eqs. (19), (24), (25) and (28) it follows that in the energy range El E <C Eh 
the photon-ALP mixing is maximal (8 ~ tt/4) and the conversion probability becomes 
energy-independent. This is the so-called strong-mixing regime. Outside this regime 
the conversion probability turns out to be energy-dependent and vanishingly small, so 
that E L and E H acquire the meaning of low-energy and high-energy energy cut-off, 
respectively. 



3.4- Photon transfer function and polarization 

As stressed in Section 1, the magnetic field crossed by the beam on its way from the 
source to us is modelled as a network of magnetic domains. For each magnetic field 
component (source, extragalactic, intracluster and Galactic) all domains are supposed 
to have the same size L equal to the coherence length and in every domain the magnetic 
field B is assumed to have the same strength but its direction is allowed to change 
randomly from one domain to another. 

As a matter of fact, the application of this approach to extragalactic magnetic 
fields neglects cosmological effects which are potentially important, since the GRB 
polarization can be measured for sources out to a redshift of about 2||. However, for 

|| Private communication from E. Costa and K. Toma. 
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the sake of clarity - and even because this effect is irrelevant outside cosmology - we 
present below the investigation of the photon-ALP mixing in the above random network 
discarding redshift-dependent complications, which will be discussed in Subsection 3.5. 

In order to accomplish our task, we take the x, y, z coordinate system as fixed once 
and for all, denoting by ip the angle between By and the y axis in a generic domain and 
treating ip as a random variable in the range < ip < 2n. 

So, what we need in the first place is the generalization of the previous result for an 
arbitrary orientation of in a single domain. Clearly, the mixing matrix M. presently 
arises from Ai^ through the similarity transformation 

M = V\iP)M^V(iP) (29) 

operated by the rotation matrix in the x-y plane, namely 



\ 

Accordingly we find [45] 



( cos ip 
sin^ 




— simp 
cosip 




0\ 




(30) 



M 



l/X 



'■yy 



V 



A a7 simp A a7 cos ip 



A a7 simp \ 

A a7 cos ip 



with 



A xx = A || sin 2 ip + Aj_ cos 2 ip , 

A xy = A yx = (Ay - A ± ) sin^ cos?/' , 



*yy 



Am cos 2 -?/' + Axsin 2 -?/' . 



(31) 

(32) 
(33) 
(34) 



Our next step consists in the evaluation of the transfer matrix in the considered 
domain. As before, our strategy is to diagonalize the mixing matrix. Dealing with Ai 
is slightly more complicated than dealing with but Eq. (29) allows to reduce the 

present problem to the one solved above [56]. It is convenient to label all quantities 
pertaining to the generic k-th domain with the index k, with the understanding that 
all ^-dependent quantities labelled by k are evaluated at the edge of the k-th domain 
closer to us. The source is located at z = and we suppose that there are iV magnetic 
domains with size L along the line of sight. Hence, the beam propagation over the k-th 
domain is described by the polarization density matrix 

Pk = T k T ] k , (35) 

where the transfer matrix T k can be rewritten as 

T k = V\i) k ) e WL W V{i) k ) , (36) 

where ip k is the random value of ip in the domain in question. The explicit form of T k 
can be represented as 



T k = e WlL T A ^ k ) + e^ L T B ^ k ) + e^ L T c ^ k ) , 



JD 2 L 



iD 3 L / 



(37) 
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with 



T B (ipk) = 



T c (ip k ) = 



V 



cos 2 Ipk 
sin ip k cos ipk 




sin 2 ip k 




sin 2 ipk cos 2 9 
sin ipk cos ipk cos 2 9 
^ sin ip k cos 9 sin 6* 

/ sin 2 9 sin 2 -0 fc 
sin 2 # sin ipk cos ^fc 



sin ipk cos -0fc \ 

, (38) 

°/ 

sin ipk cos -0fc cos 2 # sin 9 cos # sin ipk \ 

cos 2 -0fc cos 2 9 sin 6> cos 6* cos ipk 
cos ^ sin 9 cos 6* sin 2 # 

sin 2 9 sin cos ipk — sin # cos 9 sin \ 

sin 2 9 cos 2 — sin 9 cos # cos ^ 
— sin 9 cos # cos ipk cos 2 # 



(39) 



.(40) 



y — sin 9 cos 9 sin ^ 

Altogether, the whole beam propagation from the source to us is described by 

p N = T(ip N ,...,ip 1 )p T^ip N ,...,ip 1 ) , (41) 
where po is the polarization density matrix at emission and we have set 



N 



T(lpN, ■■; "Pi) = Y\T k . 



(42) 



k=i 



Finally, the 2x2 photon polarizaton density matrix (i.e. the 1-2 block of the density 
matrix for the photon- ALP system) can be expressed in terms of the Stokes parameters 
as [57] 

1 ( I + Q U-iV \ 



Pi 



U + iV I-Q 



and the degree of linear polarization Hl is defined as [58] 



n, 



(Q 2 + U 2 ) l/2 {Pll ~ P22f + (pl2 ~ P2l)' 



1/2 



(43) 



(44) 



I Pn + P22 

In the following, we will evaluate the ALP contribution to Il £ in the various physical 

situations mentioned above. 



3. 5. Effects of the cosmic expansion 

As explained in the previous Subsection 3.4, the considered strategy has to be 
slightly revised in order to account for the cosmic expansion when considering photon 
propagation in the extragalactic magnetic fields. In practice, the logic remains the 
same but distances have to be parameterized in term of the redshift z rather than 
by the coordinate along the propagation direction^. As is well known, the distance 
dl(z) travelled by a photon over an infinitesimal redshift interval dz is expressed by the 
relation 

dz 

dl (z) = : " 7Yf 2 , (45) 

H (l + z) (l + z) 2 (Q M z + l)-Q A z(z + 2) 



% In the present Subsection, z denotes from now on the redshift rather than the z-coordinate, and so 
no confusion will arise. 
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where for the standard ACDM cosmological model we have H ~ 72Kms~ 1 Mpc~ 1 
for the Hubble constant, while Qm — 0.3 and Q\ ~ 0.7 represent the average 
cosmic density of matter and dark energy, respectively, in units of the critical density 
p cr ~ 0.92 ■ 10 -29 gcm~ 3 [59]. Hence, a generic domain extending over the redshift 
interval [z a , z{\ (z a < zj,) has size 

rz b dz 

L(z a ,z b )~ 4.17/ - T7 2Gpc.(46) 

Jza (1 + z) [(1 + zf (0.3z + 1) - 0.7z (z + 2) 

After photons have been emitted by a GRB, they propagate in the intergalactic 
medium (IGM). Within the present context, it looks natural to describe the overall 
structure of the cellular configuration of the extragalactic magnetic field by a uniform 
mesh in redshift space, which can be constructed as follows. In the lack of any reliable 
information, we assume for definitiveness L Q = 1 Mpc as coherence length of the 
magnetic field B at redshift z = 0, which in our notation translates into 

L(0, zx) = 1 Mpc , (47) 

which fixes the size in redshift space of the domain closest to us through Eq. (46). 
Actually, since we are employing a uniform mesh in redshift space, Z\ sets the size of 
all magnetic domains by recursive application of Eq. (46). Hence, the n-th one extends 
from z — (n — l)zi to z — nz\ and its length can be written as L((n — l)z±, nz\). 

The absence of the Gunn-Peterson effect is usually taken as evidence that the 
IGM is ionized. The high electrical conductivity of the IGM allows us to assume that 
the electron number density n e traces the cosmic matter distribution. But since n e is 
proportional to the mass density p and the average number density of electrons n e goes 
like (1 + z) 3 , we obtain relationship 

n e (z) = (l + 5(z))n efi (l + z) 3 , (48) 

where S(z) = (p(z) — p{z))/p{z) is the mass density contrast, which quantifies the local 
deviation of the matter density from the average and can be computed starting from 
a given spectrum of density perturbations (see, e.g. [56]). Actually, two facts follow at 
once. First, the ionized nature of the IGM implies that the plasma frequency varies 
with z as 

u pl = (l + 5(z))^co plfi (l + zf 2 (49) 

owing to its dependence on n e (u> p i,o obviously corresponds to n e ,o)- Second, because 
of the resulting high conductivity of the IGM, the magnetic flux lines can be thought 
as frozen in the IGM. Therefore, flux conservation during the cosmic expansion entails 
that B scales like n^ 3 , thereby implying the magnetic field in the domain at redshift z 
is 

B = (l + 5(z)) 2 / 3 B (l + z) 2 , (50) 

where Bq denotes the average magnetic field at z — 0. 

As a rule, simulations of the mass distribution in the Universe predict that most 
of the considered domains should be under-dense with respect to the mean by a factor 
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~ 10 [60, 61]. Simulations as well as observations also show the existence of over-dense 
regions with respect to the mean by a factor 10 — 100 at low redshift [62], and Eq. (50) 
entails that in the latter regions the magnetic field should be larger by a factor ~ 5-20. 
However, the characteristic size of such over-dense regions turns out to be < 100 kpc, 
and so they fill only a small fraction of space [63]. As a consequence, the overall effect of 
these over-dense regions on photon- ALP mixing is negligible: ignoring them effectively 
amounts to take the source slight nearer and nothing else. We are therefore led for 
simplicity to discard matter overdensity effects, thereby setting 5(z) = 0. In this 
way, we also avoid committing ourselves with any specific (model-dependent) choice 
for the matter density distribution. Hence, in the following we will simply assume 
B ~ B (1 + z) 2 , where for definitiveness we will take B = 10~ 9 G as average value of 
the magnetic field at z — 0. 

Moreover, the exact value of u p \ and its scaling with the redshift [Eq. (49)] can 
also be safely neglected. The reason is that we are working within the strong-mixing 
regime, in which plasma effects are irrelevant, and so our result is independent of the 
precise value of u p \ in every magnetic domain: all what matters is that the strong-mixing 
condition is met in any domain and we have cheked that for our preferred values of the 
input parameters this is indeed the case. Finally, we remark that another consequence 
of the strong- mixing regime is that photon- ALP mixing is energy-independent, and so 
we do not have to worry about the change of the photon energy along the line of sight 
caused by the cosmic expansion. 

4. ALP effects on the polarization of gamma-ray sources 

Equipped with the results of the previous Section, we are now ready to investigate the 
implications of photon-ALP mixing on the polarization of GRBs. Needless to say, we 
also have to keep the photon-ALP conversion probability under control in order to make 
sure that the resulting dimming does not prevent the source from being detected. As 
far as the magnetic field is concerned, we will address in turn the various configurations 
crossed by the beam, namely the component inside the source, the extragalactic and the 
intracluster contributions and the Galactic component. For definiteness, we shall take 
g ai = 10" 11 GeV" 1 throughout. 

4-1. Source magnetic field 

The relevance of ALP-induced polarization in the source magnetic field for the GRB 
emission has been studied in [64]. In that paper, the authors claim a potentially large 
effect. Lacking a detailed model for the GRB environment, we follow [64] and we 
consider the source as a region where the magnetic field is homogeneous, with strength 
B ~ 10 9 G and size L ~ 10 9 cm. The typical electron number density in the source is 
estimated as n e ~ 10 10 cm -3 . Accordingly, from Eq. (27) we find the the high-energy 
cut-off (arising from QED vacuum polarization) is Eh — 2 x 10~ 3 keV, regardless of 
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m. Thus — contrary to the previous claim — in this situation ALPs do not affect the 
GRB polarization. However, this model for the source is oversimplified and in particular 
the magnetic field could have a turbulent and irregular structure. Since there are large 
uncertainties in the parameters of the GRB emitting regions, it would be premature 
to exclude definitively the possibility of any effect of ALP conversions in the source. 
However, lacking of a detailed model for the GRB emitting region, hereafter we neglect 
possible photon-ALP conversions occurring in the source. 

4-2. Extragalactic magnetic field 

Once emitted by a GRB, photons propagate in the IGM, which is modelled as explained 
in Subsection 3.5. Our results can be summarized as follows. 

In Figure 1 we show the average final linear polarization as a function of 
the photon energy E, starting with a completely unpolarized source at distance 
d = 100 Mpc. The averaging process is over an ensemble of 10 4 realizations of the 
random magnetic field network. Since in this case we are considering a relatively 
close source, we neglect the redshift dependence of the extragalactic magnetic field. In 
agreement with Eq. (26), we find that ALP effects on the polarization start to become 
relevant at E ~ 25 keV for m = 10~ 13 eV and at E ~ 2.5 keV for m = 10~ 14 eV. 
Moreover, ALP effects saturate fast as we enter the strong-mixing regime, where the 
photon-ALP conversion probability is energy- independent. A similar trend concerns the 
photon survival probability P 77 . In the following, we will focus on the strong- mixing 
regime, choosing E = 100 keV as observed reference photon energy and m = 10~ 13 eV 
as reference ALP mass. 

In Figure 2 we exhibit the average linear polarization and the photon survival 
probability P 77 as a function of the source distance d, starting from a source fully 
polarized along the x direction (p xx {0) = 1)- Again, the averaging process is over 
an ensemble of 10 4 realizations of the random magnetic field network. For our input 
values, we get A„ 7 L <C 1 in every magnetic domain, so that the conversion probability 
in a single domain dictated by Eq. (24) assumes the form P$ = (A aj L) 2 . In this 
limit, neglecting the redshift dependence in the extragalactic magnetic field, the average 
photon survival probability reduces to the simple analytic expression [45, 65] 



thereby implying on average a complete equipartition of the photon/ALP ensemble after 
passage through many domains (d >> L). As a check of our numerical simulations, in 
Figure 2 we also compare the numerical (dashed curve) and the analytical (continuous 
curve) results for P 77 (d). The nice agreement between the two prescriptions reassures 
us about the reliability of our simulations. We find that the average photon survival 
probability saturates at 2/3 for <i;>500 Mpc. Instead, the saturation of the average 
final linear polarization occurs at a larger distance d > 10 3 Mpc with asymptotic value 
U L ~ 0.79. Furthermore, an important conclusion to be drawn from Figure 2 is that 




(51) 
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Figure 1. Average final linear polarization Hl as a function of the photon energy E 
after propagation in the extragalactic magnetic field for ALP masses m = 10~ 13 eV 
(solid line) and m = 10~ 14 eV (dashed line), respectively. The emitting GRB is 
assumed to be completely unpolarized and at distance d — 100 Mpc. 



the dimming arising from photon-ALP conversion never prevents the observation of the 
gamma-ray flux. 

In Figures 3, 4, 5 and 6 we display the final probability density functions fu for 
the GRB linear polarization, assuming that all sources are at redshift z = 0.03, z = 0.3, 
z = 1 and z = 2, respectively. From our Eq. (46), these redshifts would correspond to 
distances d ~ 122.2 Mpc, d ~ 1.2 Gpc, d ~ 2.4 Gpc, d ~ 3.5 Gpc, respectively. In these 
simulations we are taking into account the redshift effects discussed in Subsection 3.5. To 
obtain the probability distributions, we have performed 10 4 simulations of the photon 
evolution within the random magnetic field configurations. For simplicity, we have 
assumed that in each simulation all sources have the same initial linear polarization, 
II = 0.0, 0.3, 0.7, 1.0, rispectively. It turns out that photon-ALP mixing smears out 
the initial GRB linear polarization. This effect increases with the source distance and 
saturates at z ~ 2. For z<>l, fu still presents a peak which is the relic of the initial 
linear polarization, but with long tails which are not present in the standard polarization 
distributions described in Section 2. In particular, the expected clustering in the linear 
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Figure 2. Average linear polarization II z, (dotted line) and photon survival probability 
P 77 evaluated numerically (dashed line) and analytically (solid line), as a function of 
the source distance d for photon- ALP mixing in the extragalactic magnetic field. The 
source is assumed to be fully polarized along the x direction. 

polarization distributions gets smeared out by photon-ALP mixing. Moreover, since 
the standard GRB linear polarization should have an initial spread and the GRBs are 
actually distributed over a variety of distances, in the presence of photon-ALP mixing 
the final distributions should to be even more smeared out and irregular than the ones 
presented in these Figures. For z^l we find rather flat probability distributions in 
the GRB polarizations, without any record of the initial linear polarization Ilo. As a 
consequence, the presence of photons- ALP mixing appears to hinder the possibility to 
extract from observational data information on the initial polarization mechanism for 
GRBs. Conversely - in the lack of any standard explanation - detection of the features 
presented above can be interpreted as a hint at the existence of very light ALPs. 

4-3. Intracluster magnetic fields 

Given that GRBs are at cosmological distances, there is a nontrivial chance that in some 
cases their line of sight crosses a cluster of galaxies. So, it is worthwhile to investigate 
what happens in this instance. 
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Figure 3. Probability density function /n for the final linear polarization II l after 
propagation in the extragalactic magnetic field, considering 10 4 GRBs at redshift 
z = 0.03, with initial linear polarization IIo = 0.0,0.3,0.7, 1.0. 



Observations have shown that the presence of magnetic fields with average strength 
B ~ 10 -6 G is a typical feature of the intracluster region. Even more remarkable is the 
fact that observations are able to yield information about the associated coherence 
length, which turns out to be L ~ 10 kpc [66]. As in the case of extragalactic 
magnetic fields, we assume a cellular structure for the intracluster fields, with domain 
size L ~ 10 kpc. As before, plasma effects are expected to be present. Specifically, 
the electron density in the intracluster medium is n e ~ 1.0 x 10~ 3 cm -3 [67], yielding a 
plasma frequency uj v \ ~ 1.2 x 10~ 12 eV. 

With these input values, the low-energy cut-off E L in Eq. (26) is E L ~ 3.6 keV for 
m < 10~ 12 eV, thereby guaranteeing that for our benchmark values E = 100 keV and 
m = 10~ 13 eV we are in the strong-mixing regime, where the photon-ALP conversion 
probability is energy- independent. 

In Figure 7 we plot the average linear polarization and the photon survival 
probability P 77 as a function of the distance D traveled inside the intracluster region, 
assuming that the source is fully polarized along the x direction. The averaging process 
is again over an ensemble of 10 4 realizations of the random magnetic field network. We 
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Figure 4. Probability density function /n for the final linear polarization II l after 
propagation in the extragalactic magnetic field, considering 10 4 GRBs at redshift 
z = 0.3, with initial linear polarization ITo = 0.0, 0.3, 0.7, 1.0. 



find that for D > 200 kpc both the polarization and the survival probability saturate at 
their limiting values, P 77 = 2/3 and lit, ~ 0.79 respectively. 

As before, the source dimming turns out to be unimportant. Thus, the result for the 
average final polarization shown in Figure 7 implies that whenever the line of sight to a 
GRB crosses a cluster of galaxies the polarization distributions shown in Figures 3, 4, 5 
and 6 are expected to be further smeared out. 

4-4- Galactic magnetic fields 

Finally, the beam from a GRB obviously crosses the Milky Way before being detected. 
Observations over the last three decades have led to a rather detailed picture of the 
Milky Way magnetic field. Perhaps, its most important feature is that it consists of two 
components, a regular and a turbulent one. We proceed to address the ALP-induced 
polarization and the photon- ALP conversion in these two magnetic field configurations. 

Regular component - Measurements of the Faraday rotation based on pulsar 
observations have shown that this component is parallel to the Galactic plane. Its 
strength varies between B ~ 2 x 10~ 6 G in the Solar neighborhood and B~4x 10~ 6 G 
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Figure 5. Probability density function /n for the final linear polarization II l after 
propagation in the extragalactic magnetic field, considering 10 4 GRBs at redshift z = 1, 
with initial linear polarization IIo = 0.0, 0.3, 0.7, 1.0. 



at 3 kpc from the center [68]. Moreover, the associated coherence length is L ~ lOkpc. 
Inside the Milky Way disk the electron density is n e ~ 1.1 x 10 -2 cm' 3 [69], resulting in 
a plasma frequency u p \ x 4.1 x 1CT 12 eV. Assuming the whole Galaxy as a single domain 
with an homogeneous magnetic field, for these input values Eq. (26) entails El ~ lOkeV 
for m < 10~ 12 eV, so that for our reference values E = 100 keV and m = 10 -13 eV we 
are in the strong-mixing regime. In Figure 8 we show the linear polarization Ul and the 
photon survival probability P 77 as a function of the photon path D in the Galaxy, taking 
as Galactic magnetic field B ~ 4x 10~ 6 G, oriented along the y direction. Both quantities 
depend on the initial photon polarization. In the left panel we have chosen the initial 
linear polarization along the y axis. In this situation, we find a final survival probability 
P 77 ~ 0.68, while the initial linear polarization II ^ = 1 remains unchanged. Instead, 
in the right panel we have started our simulation with a completely unpolarized state 
(n^ = 0). In this case, the photon- ALP conversion produces seizable effects on both the 
photon flux and the photon polarization. In particular, the survival probability after 
Galaxy crossing is P 77 ~ 0.84 and the final linear polarization il^ ~ 0.19. A detailed 
study of this effect would require a better modelling of the Galactic magnetic field, as 
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Figure 6. Probability density function /n for the final linear polarization II l after 
propagation in the extragalactic magnetic field, considering 10 4 GRBs at redshift z = 2, 
with initial linear polarization IIo = 0.0, 0.3, 0.7, 1.0. 



in the study of photon- ALP conversion in the Milky Way performed in [70]. However, 
we leave further refinements of these estimates for future work. 

Turbulent component - Over much smaller scales, the dominant Galactic magnetic 
field appears to be stochastic, with a Kolmogorov spectrum a = 5/3 [71]. In practice, 
this component can be described by a cellular structure, with strength B ~ 1 x 10~ 6 G 
and domain size L ~ 10~ 2 pc. It is straightforward to realize that in this case the 
oscillation length l osc = 27r/A osc is much larger than the domain size L, so that the 
photon- ALP conversion is vanishingly small. 

5. Conclusions 

Measuring the polarization of prompt GRB emission in the keV-MeV range represents 
one of the main challenges for high-energy astronomy of the next decade. Various 
polarimetric missions are currently being developed, which are expected to collect an 
all-sky rich sample of GRBs so as to allow for a meaningful statistical analysis of their 
polarization properties. As recently realized, important conclusions concerning the GRB 
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Figure 7. Average linear polarization ITz, (dotted line) and photon survival probability 
P 77 (dashed line) as a function of the distance D traveled inside the intracluster region. 
The source is assumed to be fully polarized along the x direction. 



emission models are expected to be drawn from these studies. 

We have shown that the existence of ALPs with parameters lying in experimentally 
allowed ranges drastically modifies the GRB polarization pattern. More specifically, 
cosmic magnetic fields of extragalactic, intracluster and Galactic origin along the line 
of sight to a GRB act as catalysts for significant photon- ALP mixing. In particular, 
due to the random structure of the extragalactic and intracluster magnetic fields, the 
amount of photon- ALP mixing strongly depends on the orientation of the line of sight. 
Therefore, starting with a given source polarization a broad statistical distribution is 
expected to be detected when observing GRBs from different directions in the sky. 
This effect should be superimposed on the original distribution of the GRB polarization 
as predicted by standard emission models. As a result, the scatter in the distribution 
arising from the photon- ALP mixing hinders the possibility to extract information on the 
initial polarization for the observed GRBs. Alternatively, in the lack of any standard 
explanation, the observation of these peculiar broad distributions would hint at the 
existence of very light ALPs. We have restricted throughout our attention to the linear 
polarization because in all cases that we have analyzed the circular polarization turns out 
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Figure 8. Linear polarization IIx (solid curve) and survival probability P 11 (dashed 
curve) for the photon- ALP mixing in the regular component of the Galactic magnetic 
field as a function of the photon path D in the Galaxy. We have assumed photons to be 
fully polarized along the y direction (left panel) or completely unpolarized (right panel) 
before entering the Galaxy. We have taken as Galactic magnetic field B — Ax 10~ 6 G, 
oriented along the y axis. 



to be vanishingly small. This circumstance is at variance with the case of chameleontic 
ALPs [35], and so it can be used to further investigate the nature of ALPs. 

In the last few years, different hints and possible signatures of very light ALPs 
have been proposed in connection with astrophysical observations [54, 55, 65, 70, 72, 
73, 74, 75, 76, 77, 78]. In general, the presence of ALPs can play a role for values of 
the photon-ALP coupling constant g ai < 10 -11 GeV -1 , namely one order of magnitude 
lower than the current experimental limit set by the CAST experiment. This indicates 
the high potential of astrophysical observations to probe ALPs in a region of the their 
parameter space otherwise unreachable. It remains to see whether the above hints at 
ALPs will ultimately converge in a robust signature. Every astrophysical argument has 
its own systematic uncertainties and its own recognized or unrecognized loop-holes, 
so that in order to corner the ALPs it is important to use as many independent 
channels and as many different approaches as possible. In this sense, we believe 
that the present study of ALP effects on the GRB polarization should be regarded 
as an additional opportunity to uncover the existence of very light ALPs through 
astrophysical observations. Remarkably enough, no dedicated experiment is required 
and our predictions can be easily tested in the near future with the planned polarimetric 
missions for GRBs. 
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